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Abstract To reproduce the weak magnetic field on the polar caps of the Sun observed 
during the declining phase of cycle 23 poses a challenge to surface flux transport models 
since this cycle has not been particularly weak. We use a well-calibrated model to 
evaluate the parameter changes required to obtain simulated polar fields and open flux 
that are consistent with the observations. We find that the low polar field of cycle 23 
could be reproduced by an increase of the meridional flow by 55% in the last cycle. 
Alternatively, a decrease of the mean tilt angle of sunspot groups by 28% would also 
lead to a similarly low polar field, but cause a delay of the polar field reversals by 1.5 
years in comparison to the observations. 
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1 Introduction 

The current solar minimum shows some peculiar features. There were 265 and 261 spot- 
less days in the years 2008 and 2009, respectively. Only the years 1878, 1901 and 1913 
showed more spotless days in the daily records since 18493- The polar field is about 40% 
weaker than during theprevious three minima for which routine magnet ograph mea- 
surem ents are availablqj. Polar coronal-hole areas shrank by about 20% ( Wang et al.l 



l2009h and the heliospheric open flux reached the lowest values since systematic mea- 
surements started in 1967 (|Wang et al.ll2009h . Observations from Ulysses showed 



that the high-latitude solar wind is 17% less dense and cooler than during the previous 
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solar minimum i McComas et al.ll2008h . The unusual properties of the polar corona, 
the open flux, and the solar wind are probably all closely related with the weakness of 
polar field. 

This situation poses an interesting challenge to surface flux transport models (SFTM), 
which consider the passive transport and evolution of the radial component of the 
magnetic field on the solar surface under the effect of differential rotation, merid- 
ional flow and diffu s ion JWang et a.1.11 1989): Ivan Balleeooiien et aljEiii ; ISchriiverlEoOll : 
iMackav et ai1l2002l :l Baumann et al.ll2004l . 120061 ). Such models have been shown to be 
consistent with the observed evolution of the magnetic field distribution at the sur- 
face, including the reversals of the polar fields. The models take their input from the 
observed flux emergence (e.g., in terms of sunspot data or magnetograms) , so that 
their results are closely tied to the activity level of the solar cycles. Since cycle 23 was 
not particularly weak in comparison to its predecessors, the polar fields generated by 
the flux transport models tend to be much stronger than observed. In this paper, we 
consider which changes of the parameters or input data would be required for cycle 23 
in order to bring the models into agreement with the actual data. 



2 Solar surface flux transport model 



In a preceding study ( Cameron et al.| [2010). we have carried out flux transport simu- 



lations including the observed cycle-to-cycle variation of sunspot group tilt angles to 
model the solar surface field and open flux from 1913 to 1986, based upon the RGO and 
SOON sunspot area data sets0. We obtained a reasonable agreement of the time evolu- 
tion of the open flux and of the reversal times of the polar fields with empirical results. 
Unfortunately, this model cannot easily be extended to cover the last cycles because 
systematic observations of tilt angles are not available and the definition of sunspot 
groups is inconsistent between the RGO and SOON data sets. We have therefore used 
the monthly international sunspot number since 1976 as the basic input quantity de- 
scribing the eme rgence of new flux (c f. Fig. [TJ. This method has been validated by 
compa rison with lCameron et al ] l|2O10l ) from 1913 to 1986 (for details, see I Jiang et al.l 
1201 lbl ). 

We take the number of new bipolar magnetic regions (sunspot groups) during a 
given month to be proportional to the corresponding sunspot number. The SFTM 
then requires the prescription of latitude, longitude, area, and tilt angle for each new 
bipolar magnetic region. This is done on the basis of empiri cal relations with cycle 
strength together with using random numbers (for details, see ljiang et al-lfeoilah . As 
an example, Figure [5] shows the time-latitude distribution of emerging bipolar regions 
as model input for the three cycles considered here in comparison with the actually 
observed butterfly diagram. 

The area of bipol ar regions is chose n randomly according to the observed size 
distribution function I Jiang et ah 2011a!'). with a higher probability for large groups 



during cycle maximum phases (lHathawavl I2010T ). Since preferred longitudes of flux 
emergence (activity nests) strongly affect the equatorial dipole component and thus the 
open flux during solar maxima, we draw the emergence longitudes from a combination 
of random longitudes with a set of longitudes clustering on two values 180 degrees 
apart (for details, see Ijiang et al.|[2011all2011rj ). On the basis of Mount Wilson and 



3 http:/ /solarscience. msfc.nasa.gov/greenwch.shtml 
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Fig. 1 Monthly international sunspot number since 1976. 



Kodaikanal sunspot data. iDasi-Esrmig et alJ ( 2010h found that the average tilt angle 
of sunspot groups is negatively correlated with the strength of the cycle. The cycle- 
dependent factor T n in the fit relation a n (A) = T n y/\ (a n : average tilt angle for cycle 
n; A: latitude) is determined by using this anti-correlation ((Carneror^ i etaL||20^j). The 
values of T„ for n = 21, 22, and 23 are 1.21, 1.21, and 1.32 (jjiang et al.ll2011bFT 

The profile for the differential rotation used in our SFTM simul ations is {2(A) = 
13.38 — 2.30 sin 2 A — 1.62 sin 4 A (in degrees per day; Snodgrass|[l983h . The profile for 
the meridional flow is (cf. Ivan Ballegooiien et al.lfl998T ) 



v(X) 



llsin(2.4A) ms~ 




whcrc|A| < 75° 
otherwise, 



(1) 



which is shown by the solid line in Figure El The supergran ular turbulent diffusivity is 
chosen as 250km 2 s -1 (following ISchriiver and Zwaanll2000l . Table 6.2). 

The i nitial magnetic field d istribution at the start of the simulation is given by 
Eq. 10 of I Cameron et al. | <|2010h (see also van Ballegooijen et al. 1998). The amplitude 
of the initial field has been set such that the line-of-sight field averaged over a 35° 
wide polar cap agrees with the corresponding WS Q data for 1976 (see F ig. 5) . A decay 
term due to radial diffusion of the magnetic field ( Baumann et al. 120061 ') has not been 
considered her e since the observed an ti-correlation of sunspot group tilt angle and 
cycle strength dDasi-Espuig et al. 2010h removes the necessity for this term in short- 
term studies dCameron et al. 2010T ). These choices of the model parameters and the 
flux input constitute our 'reference model'. 



3 Results of the reference model 

Figure [3] shows the time evolution of the simulated radial component of the polar 
fields, averaged over 15°-wide polar caps. Since the location, size and emergence time 
of the bipolar magnetic regions involve random numbers, the results from different sets 
of random numbers vary somewhat. The error bars in Figure |3] denote the standard 
deviation for 20 runs with different random numbers. In order to compare with the line- 
of-sight polar field above 55° latitude observed at Wilcox Solar Observatory (WSO), 
we use the same definition of the polar field as WSO. The result is shown in Figure[5] 
The WSO polar field is multiplied by a constant factor 1.3 to correct for magnetograph 
saturation ( Svalgaard et al.lll978T ). There is good agreement between the simulation and 
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Fig. 2 Time-latitude distribution (butterfly diagram) of the bipolar magnetic regions used as 
input for the SFTM (lower panel) in comparison to the actual sunspot butterfly diagram from 
the SOON data (upper panel). 




Latitude 

Fig. 3 Profile of meridion al flow adopted in our reference model (solid curve) and the profile 
used by Wang ct al. (2009, dashed curve). 



the measurement for both the amplitude and reversal times of the polar field between 
1976 and 2002. However, during the end of cycle 23, the simulation gives polar fields 
that are too strong by a factor of about 2. 

In addition to the polar field, we also use the heliospheric open flux to compare 
observation and SFTM simulation. To this end, we extrapolate the simulated photo- 
spheric fi*U_out_t£_the_s^urce_sm source surface (CSSS) 
model dZhao and Hoeksemalll995blll995al ; |jiang et al.ll2010l ) with the model parameters 



given by Cameron et al. ( 2010h . Figure [6] shows the comparison betwe en the simulation 
results and the OMNI spacecraft data (with a kinematic correction, cf . lLockwood et aD 
l2009al ). 

The simulated open flux is consistent with the measurements in both maximum and 
minimum phases before 2005. During the minimum of cycle 23, however, the simulated 
result is much higher than the measurement. This is directly related to the too strong 
polar field since the open flux during activity minima is dominated by the low-order 
axial multipoles. 
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Fig. 4 Simulated evolution of radial component of solar north (full line) and south (dashed 
line) polar field, averaged over the 15°-wide polar caps. The error bars denote the standard 
deviation for 20 runs with different sets of random numbers. 




1990 1995 



Fig. 5 Comparison between the simulated (red curves) and observed (WSO data, multiplied 
by a factor 1.3, black curves) linc-of-sight component of solar polar field, averaged over 35°- 
wide polar caps. Full lines correspond to the north pole, dashed lines to the south pole. 




Fig. 6 Evolution of heliospheric open flux since 1976. Shown is the result of the reference 
SFTM (red curve) in c omparison to the OMNI spacecraft data (with kinematic correction, 
Lockwood ct al . 2009a. black curve) and va lues derived from geomagnetic variations (with 
kinematic correction. lLockwood et al.ll2009bl , grey curve). 
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4 Possible cause of the weak polar field in cycle 23 

What could have been different in cycle 23? Which changes in the SFTM parameters 
and/or input properties could possibly lead to the low polar field? The turbulent dif- 
fusivity models the random walk of magnetic elements due to supergranular motions. 
Sinc e these motions are observ ed to only weakly vary over the solar cycle (for a review 
see lRieutord and Rincon1l20ld . section 4.6.5), we have assumed that the diffusivity is 
constant in time. Likewise, the time-latitude distribution of emerging flux during cycle 
23 is very similar to that of the preceding cycles (cf. Fig. [5| • The same is true for the 
mean sunspot group areas as derived from the SOON data . 

We consider three possi b ilities for potential parameter variations, a) sunspot num- 
ber: IWilson and Hathawavl l|2005h found that the international sunspot number has 



been systematically overestimated since 1981. Lower values of the sunspot number 
would entail le ss flux emergence i n our model and thus a lower polar field, b) tilt 
angle: although ISchrii ver and Liu (|2008f ) found no indication for a systematic change 



in the mean tilt angle by analyzing selected MDI magnetograms, one has to keep in 
mind that the det ermined tilt angles are sensitive to the location and evolution phase 
( Yang et all 120091 ) of the sunspot group considered, to the observer's definition, and 



to the instrument. Furthermore, the tilt angles show a very large scatter, possibly due 
to th e convective buffeting of rising flux tubes in the upper layers of the convection 



zone ( Longcope and Choudhurill2002h . c) meridional flow: significant var i ability of the 



meridional flow has been detected during cycle 23 dChou and Pail l200ll : iGizonI 120041 : 
lHathawav and Rightmirell2010l : iGonzalez Hernandez et af I l2010h . We considered vari- 



ations of the meridional flow, of the mean sunspot group tilt angle, and of sunspot 
number during cycle 23 as potential causes of the low polar field during that cycle and 
used our SFTM to determine which changes of these properties in cycle 23 would be 
required to reproduce the observations. 

Figures [7] and [8] show the evolution of the polar field and the heliospheric open 
flux, respectively, that result from such changes (varying one parameter at a time): a) 
a 40% reduction of the sunspot number (green curves), b) a 28% decrease of mean tilt 
angle (blue curves), or c) a 55% increase of the meridional flow (red curves), always 
with respect to the values for cycle 23 in the reference model. In case a), the open flux 
during the maximum phase of cycle 23 is much weaker than the observed value since 
the decreased amount of flux emergence leads to a strongly reduced equatorial dipole 
moment during the maximum phase. This excludes an overestimate of the sunspot 
number as sole cause of the low polar field. For case b), the reduced tilt angles lead 
to a delay of the polar field reversal by about 1.5 yrs, which is in disagreement with 
the observations. Smaller tilt angles correspond to less meridional magnetic flux of the 
bipolar magnetic regions, so that it takes longer (more bipolar regions) to reverse the 
polar field of the preceding cycle. 

For the case c) with increased meridional flow, both polar field reversal time and am- 
plitude match the observation (cf. Figure [7} . A stronger meridional flow, especially an 
increased velocity gradient near the equator, effectively separates the two hemispheres, 
so that the transport of net flux to the poles is reduced. Although the contribution 
of each magnetic region to the polar field is smaller, the transport velocity is higher, 
thus minimizing the effect on the polar field reversal time. The increased meridional 
flow also brings the simulated heliospheric open flux into reasonable agreement with 
the observation during both the maximum and the minimum ph ases of cycle 23 (cf . 
Figure [8]). This result is consistent with the SFTM simulations of I Wang et all (2009) 
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Fig. 7 Evolution of line-of-sight polar field with modified parameters during cycle 23: 40% 
decrease of sunspot number (green curves), 28% decrease of the sunspot group tilt angle (blue 
curves, which nearly lie upon the green curves), 55% increase of meridional flow (red curves). 
The three simulated cases are identical before the start of cycle 23. The WSO observations are 
given by the black curves. In all cases, full lines correspond to the north pole, dashed lines to 
the south pole. 




Fig. 8 Evolution of the heliospheric open flux with modified parameters during cycle 23. The 
colored lines refer to the same cases as in Figure [7] The black and grey lines give the OMNI 
data and values derived from geomagnetic variations, respectively. 



and ISchriiver and Liul ( 20081 ) . IWang et "al", | (|2009h found that in their model already a 
15% increase of the meridional flow, much smaller than our value of 55%, is sufficient 
to reproduce the observed weak polar field. This is probably due to the extremely steep 
velocity gradient near the equator in the meridional flow profile adopted by these au- 
thors (cf. Figure [3]). However, such a prof ile does not seem to be very realistic in the 
light of helioseismic measurements (e.g.. lGizon and Birchll2005n . 



5 Conclusion 



Our simulations indicate that quite substantial changes in the SFTM parameters are 
required to obtain the observed weak polar field during cycle 23, at least if the effect is 
ascribed to a change in only one parameter, such as a 28% decrease of the mean sunspot 
group tilt angle or a 55% increase of the meridional flow (from 11 ms _1 to 17 ms -1 ). 
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Note that a combination of changes in more than one parameter would of course 
require smaller variations to bring the SFTM simulatio ns into agreement with the 
observations. For instance, the results of lHathawavl 1 20ld ) indicate that the meridional 
flow speed during the decay phase of cycle 23 was nearly 20% higher than during the 
corresponding phase in cycle 22. Furthermore, the meridional flow profile used in the 
flux tran sport simulation probably unde restimates the poleward velocities beyond 60° 
latitude ( Hathaway and Rightmirg|20f ll ). On the other hand, polar countercells of the 
meridional flow could redistribute flux from the polar ca p to lower latitudes and thus 
also contribute to weaker polar fields I Jiang et al.l 12009). When additional empirical 
data (e.g., concerning the tilt angles) become available, it will be possible to further 
constrain the simulation parameters and hopefully achieve a better understanding of 
this rather unusual solar minimum. 
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